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We introduce a switching mechanism in the asymptotic occupations of quantum states induced 
by the combined effects of a periodic driving and a weak coupling to a heat bath, ft exploits one of 
the ubiquitous avoided crossings in driven systems and works even if both involved Floquet states 
have small occupations, ft is independent of the initial state and the duration of the driving. As 
a specific example of this general switching mechanism we show how an asymmetric double well 
potential can be switched between the lower and the upper well by a periodic driving that is much 
weaker than the asymmetry. 
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I. INTRODUCTION 

The interplay between a coherent periodic driving force 
and the incoherent damping of a thermal environment 
enriches the dynamics of a quantum system and opens 
new potential applications In addition to controlling 
the transient dynamics, e.g. with respect to tunneling [2|, 
H, H[, a control of the asymptotic state is desirable. The 
ability to design a system's probability distribution, e.g. 
to switch between two macroscopically distinguishable 
states, in the presence of a thermal environment is a key 
to quantum control techniques. 

The paradigmatic model for switching is a double 
well potential which is experimentally realized in super- 
conducting quantum interference devices (SQUIDs) 
atom-optical potentials @, 0], spin tunneling in con- 
densed matter [8|] or in the transfer of protons along 
chemical bonds Q. In some cases the model can be 
restricted to a two-level system. Different approaches 
for switching by a population inversion in driven two- 
level systems have been proposed, e.g. induced by 
symmetry-breaking [nj, structured environments [Til 
strong nonequilibrium noise fl2| . or strong driving [l3| . 
However, the restriction to a two- or a three-level system 
limits the possible switching mechanisms. 

Time-periodic quantum systems are best described by 
Floquet states, which are solutions of the Schrodingcr 
equation without the coupling to the environment. When 
coupling the time-periodic system weakly to a thermal 
bath, all Floquet states are asymptotically populated 
with occupation probabilities, which can be determined 
within a Floquet-Markov approach [ll [ll [H [l7l. [T^|. 
These occupations are quite different from the canonical 
distribution of Boltzmann weights in undriven systems 
and so far lack an intuitive understanding. 

In this paper we demonstrate a dramatic property of 
time-periodically driven quantum systems weakly cou- 
pled to the environment: The asymptotic state can be 
switched to an almost orthogonal state by a small pa- 
rameter variation. This is in stark contrast to time- 
independent systems, where the asymptotic occupations 
are determined by Boltzmann weights and vary slowly 
with a parameter. The proposed switching mechanism 



exploits one of the ubiquitous avoided crossings in driven 
systems and works even if both involved Floquet states 
have small occupations. As a specific example of this 
general switching mechanism we show for an asymmetric 
double well potential, see Fig. [TJ that a weak periodic 
driving switches the cycle-averaged asymptotic probabil- 
ity density from the ground state of the undriven system 
in the left well to the right well. Note, that the periodic 
driving is much weaker than the asymmetry, see Fig.[lja), 
and therefore this switching is unrelated to previous stud- 
ies on hysteretic switching in a driven dissipative double 
well [l9l , |20| . We explain the switching mechanism by an 
effective rate equation, which combines the effects of the 
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FIG. 1: a) Asymmetric double well potential and its eigenen- 
ergies without driving, A — (solid line), and the almost in- 
distinguishable variation of the potential for a small driving 
amplitude Aq « 0.008 (dashed and dotted line), b) Asymp- 
totic probability density p(x) for A = and a small tem- 
perature 1//3 = 1/100, with almost all probability in the left 
well, c) Cycle-averaged asymptotic probability density p(x) 
according to Eq. [5] for A = Aq, with more than 99 % of proba- 
bility in the right well, demonstrating a weak driving induced 
switching to a macroscopically different state. See Fig. [2] for 
parameters. 
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coherent driving at an avoided crossing of two Floquet 
states with the incoherent bath coupling. 

The paper is organized as follows: In Sec. Hi] our model 
for the periodically driven, dissipative double well is in- 
troduced. The switching mechanism is investigated in 
Sec. HID We finally conclude and discuss advantages of 
the switching mechanism in Sec. IIVI 



II. THE MODEL SYSTEM 

As an example we study a particle in an asymmetric 
double well potential in the quantum regime, where the 
ground state is in the left well and the first excited state is 
in the right well, see Fig.[lja). It is driven by an additive 
time-periodic force, leading to the system Hamiltonian 
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where p is the asymmetry parameter of the double well 
potential and A and f2 are the driving amplitude and 
frequency, respectively. We introduce the dimensionless 
quantities x = x/xq, H s = H s /Vq, t = t ■ Vq/%, £1 = 
Q, ■ h/Vo, and H e g = h/ (y/mVoXo) . In the following we 
omit the tilde and then the dimensionless Hamiltonian 
reads 



fr 2 r) 2 



+ x 4 - x 2 + x (p + AcosQt) . (2) 



The Schrodinger equation of a periodically driven quan- 
tum system has according to the Floquet theorem solu- 



tions of the form ipi(t) 



Ui{i), with ui(t + T) 



Ui{t) and T = 2-k /Vt the period of the driving. The time- 
periodic parts Ui (t) of the Floquet states form a complete 
orthonormal set at all times. The quasienergies Ei can be 
chosen to lie in the interval [0, il). 

The coupling to a heat bath is modeled in a standard 
way by a Hamiltonian [2l| 



H(t) = H s (t) + H b + H ; 
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The bath Hamiltonian Hi, = ^2 n , 

scribes an ensemble of nonintcracting harmonic oscil- 
lators coupled via H s b = x c n x n to the system. 
The properties of the system-bath coupling are given 
in terms of the spectral density of the bath J(u>) := 

f J2n m„^„ ( w ~ u n) - 5(uj + u) n )]. In the continuum 
limit the spectral density is assumed to be a smooth 
function which is linear for an Ohmic bath. An ex- 
ponential cutoff beyond the spectral mode oj c leads to 
J{lo) = e - !"'/^, where r] is proportional to the clas- 
sical damping coefficient. 

In the presence of the heat bath the state of the sys- 
tem is described by the reduced density operator p{t). 
Its equation of motion for time-periodic quantum sys- 
tems has been derived within the Floquet-Markov ap- 
proach 0, [l5|, [l6|, Qj], HH : Herein the Floquet formalism 



ensures a non-pcrturbative treatment of the driven sys- 
tems coherent dynamics. The coupling to the heat bath 
is treated perturbatively, which is valid in the limit of 
weak coupling between the driven system and the bath. 
This approximation requires a rapid decay of bath corre- 
lations compared to the typical relaxation time of the 
system and we further require CI <^ lu c . In the fol- 
lowing we restrict the discussion to the limit of large 
times, larger than the relaxation time. In this limit 
the density matrix in the basis of the periodic parts 
Ui(t) of the Floquet states is approximated as time- 
independent [3, El- Note, that the corresponding den- 
sity operator, Y^i j \ u i(t)) Pij{ u j{t)\i 1S time-periodic be- 
cause of the inherent time-dependence of the Ui(t). The 
matrix elements pij obey the rate equation 

i(e i -e j )p ij = (4) 

- ^j/'/, />'./,://. + PiiR*jkdk - Pkl [Rlj-M + Rti;lj) }■ 
k,l 

The complex rates Rij-u = ^Y Jm Xij{m)x* kl {m)g{ei - 
Ek — mil) describe bath-induced transitions between the 
Floquet states, the Xij(m) are the Fourier coefficients 
of the time-periodic matrix elements (v,i(t)\x\uj(t)) , and 
g{ijj) is the correlation function of the bath coupling op- 
erator. The latter is given by g(uS) = np(oS) J(u>)/n with 
the spectral density J(lo) and the thermal occupation 
number np{ui) of the boson bath with temperature 1//3. 
In numerical studies of the rate equation ^ one has to 
use a finite basis of Floquet states. The validity of this 
approximation is discussed in Ref . [l8j . 



III. THE SWITCHING MECHANISM 

We will demonstrate the switching process by studying 
the asymptotic spatial probability density averaged over 
one period of the driving 
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It can be expressed in terms of the solutions of Eq. ((4]) 

by P( x ) = y.,., P:,T .[,! i]i u*(x,t)ui(x,t). Figure [T^b) 
shows that for the undriven double well, A = 0, in ther- 
modynamic equilibrium at low temperatures almost all 
probability is in the left well. This reflects the dominant 
occupation of the ground state. Figure [TJc) shows that 
for a small driving amplitude, Aq w 0.008, the probability 
density is almost completely transferred to the right well. 
Note, that the driving amplitude is so small, that at all 
times the right well is energetically higher than the left 
well. This example demonstrates that a weak periodic 
driving not only alters the static Boltzmann occupation 
probabilities [l5|, but can switch to an almost or- 
thogonal and macroscopically different asymptotic state 
of the system. 

We get a first insight into this dramatic phenomenon 
from Figs. H£a) and (b), where one can see that under 
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FIG. 2: a) Quasienergy spectrum for the 17 lowest Floquet 
states vs. driving strength A and magnification of the avoided 
crossing at A — Aq (solid lines) with A = |e2 (^4o) — £7 (^4o) | ~ 
1.82 -10 -6 , the quasienergies corresponding to diabatic states 
2 and 7 (dashed lines), and eigenenergies of the undriven po- 
tential (inset), b) Stationary occupations pa in the diabatic 
basis (solid lines) and approximation based on effective rate 
R ac , Eqs. (j6]) and Q (dotted lines), c) same as b) with loga- 
rithmic axis for pa. The parameters are p — 0.03, h e g = 0.04, 
fi = ft cf f/0.768, /3 = 100, r/ = 10~ 4 and uj c = 100. 



The diagonal density matrix elements pa in the dia- 
batic basis are shown in Fig. H^b) an d (c) ■ One observes 
that poo, which corresponds to being in the ground state 
of the undriven system, drops from close to one to al- 
most zero for A = Aq. In contrast, the probability p\\ 
increases almost to one, which corresponds to the first 
excited state being dominantly populated. The tiny oc- 
cupations P22 and P77, i.e. the probabilities to be in one of 
the states of the avoided crossing, become equal. These 
observations for pa are consistent with the spatial prob- 
ability density observed in Fig. Wlc) and can indeed be 
exploited for a switching between the wells: Tuning the 
driving amplitude from outside the avoided crossing into 
its center is accompanied by a probability transfer from 
the former ground state in the left well to the first excited 
state localized in the right well. 

While the equality P22 — P77 at the center of an avoided 
crossing of states 2 and 7 is quite plausible, the main 
question is still unanswered: How can states and 1, 
which are not involved in the avoided crossing, inter- 
change their probability? 



A. Effective rate equations 

We will answer the above question by using an effective 
approximate rate system introduced in Ref. (l8| . which 
is derived from Eq. 0, 



— — Pii ^ Rik + ^ PkkRki 



(6) 



for the diagonal elements pa in the diabatic basis with 
an additional rate 



the variation of the driving amplitude A the quasienergy 
spectrum shows around A = Aq an isolated avoided 
crossing of the states 2 and 7 originating from the second 
and the 7th excited state of the undriven system. We em- 
phasize, that both the ground state, which is dominantly 
populated at A = 0, and the first excited state, which 
will turn out to be dominantly populated at A = Aq, are 
not involved in this avoided crossing. 

An intuitive understanding of the switching from the 
rate equation seems impossible: Tuning through an 
avoided crossing of the two Floquet states 2 and 7 they 
exchange their character and thus drastically affect in 
Eq. ([3]) a large number of rates Rij-ki, where one of the 
four indices is 2 or 7. In order to visualize the changes of 
the density operator due to this avoided crossing it is con- 
venient to express this operator in a basis that does not 
significantly change in the neighborhood of the avoided 
crossing. In the subspace of the Floquet states of the 
avoided crossing we use the diabatic states 2 and 7, which 
would correspond to an exact crossing. Due to the weak 
driving amplitude A <C p they are nearly identical to the 
cigenstatcs of the undriven system (Fig.^b), inset). We 
will denote quantities in this diabatic basis by a bar. 



R a 
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(7) 



replacing the rates R27, R72 in Eq. ((6|) due to the single 
isolated avoided crossing of diabatic states 2 and 7. Be- 
fore we make use of these equations, we make a number 
of remarks: The rates Rik = Rik\ik are expressed in the 
diabatic basis. The rate V = T2 + Tj + R22 + R77 — 2R22:77 
with Ti = J2k=£i R< ik (* = ^' 7) describes the transitions 
from the states of the avoided crossing to all other states. 
It is proportional to rj with a factor that is specific to 
an individual avoided crossing. The rate i? ac depends 
on the minimal splitting A of the avoided crossing and 
the dimcnsionless distance d :— (e 7 — £2)/ A from the 
avoided crossing. The main assumptions used in the 
derivation (l8j is that all quasienergy splittings ey, apart 
from the isolated avoided crossing of interest, are much 
larger than the rates Rij-ki- This is fulfilled for a suffi- 
ciently weak coupling to the heat bath and allows for ne- 
glecting almost all off-diagonal density matrix elements. 
The only non-negligible off-diagonal elements are P27 and 
P72, which are decoupled from Eq. ^ and proportional 
to P22 — P77- The dotted lines in Figs. [Hb) demonstrate 
this approximation. 
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The main advantage of the effective rate system in the 
diabatic basis, Eq. ([5]), is, that tuning the distance d from 
the avoided crossing affects exclusively the rate i? ac . In 
the center of the avoided crossing, d = 0, and for a small 
enough coupling to the heat bath, T <C A, it is much 
larger than all other rates. This leads directly to almost 
equal occupations of the diabatic states involved in the 
avoided crossing, p~22 — P77- We explain the dominant 
occupation of state 1 as the combined result of the fol- 
lowing facts: (i) The overall stationary probability flux 
between any two states is in general nonzero, as detailed 
balance is broken by the periodic driving, (ii) The rates 
between neighboring states localized in the same well are 
much larger than other intra-well rates as well as inter- 
well rates (due to the small spatial overlap between the 
states of different wells). Therefore, among the states 
confined to the same well detailed balance approximately 
holds true, e.g. between state and 2 or state 1 and 3. 
(iii) i? ac is the dominant rate and induces occupation 
equality of states 2 and 7. Figure OJc) shows the result- 
ing depopulation of state 2 towards state 7. (iv) Due to 
the approximate detailed balance among the states in the 
left well the relative occupation of state and 2 remains 
constant and therefore poo drops down together with p~22 ■ 
(v) The states in the right well equilibrate as before but 
with increased weights due to probability conservation. 
This explains the switching process observed in Fig. [TJ 

An additional surprising phenomenon is observed in 
Fig. EJb) . The impact of the avoided crossing on the oc- 
cupations occurs within a significantly broader range of 
the driving amplitude A compared to the width of the 
avoided crossing. For the parameters of Fig. [2] the full 
width at half maximum of P r {A) = J Q dxp(x), the prob- 
ability to be in the right potential well (Fig. [3](a)), is a 
factor of 30 larger than the width of the avoided crossing. 
According to Eqs. © and 0, the occupations change, 
if the magnitude of R ac is larger than or comparable to 
other significant rates in Eq. ([5]). Since these rates vary 
over many orders of magnitude, this criterion may be 
fulfilled even beyond the avoided crossing, \d\ > 1, qual- 
itatively explaining the enlarged width of P r (A). 
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FIG. 3: a) Total probability in the right well P r vs. driving 
amplitude A. b), c) Peak height P r (Ao) vs. effective coupling 
strength F/A and temperature Diamonds indicate the 

parameters of Fig. [2] The dashed line in c) gives the prob- 
ability in the right well without driving, A = 0. The dotted 
line at 1//3 = Ei — Eo indicates the transition between the 
high and the low-temperature regimes. 



pling r < A we have a high switching efficiency and one 
can show that it is independent of T in the limit Y — > 0. 

Figure[3](c) shows the influence of the temperature 1//3, 
which can be related to the level spacing E\ — E of the 
undriven system. At high temperatures, 3> E\ — Eq, 
the Floquet states are almost equally occupied resulting 
in P r (A ) ~ 0.5. For temperatures 1/(3 < E\ — E the 
probability in the right well becomes dominant (while, 
of course, it vanishes in the undriven case ^4 = 0). For 
even lower temperatures, however, P r (Ao) drops to zero. 
Here, the occupation equality p~22 — 977 is rendered by an 
increase of P77 towards P22 , hi contrast to the decrease of 
P22 towards P77 in Fig. [2jc) . Together with p~22 also poo 
remains constant and therefore switching does not take 
place. The origin of this low-temperature dependence 
remains open. 



B. Parameter dependence 

What are the optimal parameters for this switching 
effect? A maximal switching efficiency is achieved by a 
high value of the probability in the right well P r (A) = 
L°° dxp(x). This quantity is shown in Fig. Oa). Fig- 
ure [3Jb) demonstrates that if the coupling to the heat 
bath is larger than the minimal splitting of the avoided 
crossing, T > I00A, almost no probability is switched 
to the right well. (Note, that even for the largest val- 
ues of T in Fig. [3f b) the assumption of weak coupling of 
the Floquet-Markov approach is still fulfilled.) This is 
due to the fact that in this limit i? ac becomes negligible 
compared to the other rates and thus the influence of the 
avoided crossing vanishes (l8j . In contrast, for small cou- 



C. Minimal example 

A minimal example, where one of the partners of the 
avoided crossing is the ground state in the left well, is 
shown in Fig. 2] The above discussion then simplifies 
since step (iv) is eliminated and essentially just three 
states are involved. The parameter dependence of P r (Ao) 
on the effective coupling strength Y / A is unchanged and 
the switching mechanism is maintained even for low tem- 
peratures, see Fig. 

For a clear presentation we have above chosen exam- 
ples in the limit of a small driving amplitude, where the 
Floquet states are not very different from the eigenstates 
of the undriven system. In this case an avoided crossing 
requires near-resonant driving, Ei ~ E% ~ 317 in Fig. [5] 
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FIG. 4: a) Quasienergy spectrum for the ten lowest Floquet 
states vs. driving strength A and magnification of the avoided 
crossing at A — Aq (solid lines) with A = |eo(Ao) — £3(^0) | ~ 
3.57 -lO" 6 , the quasienergies corresponding to diabatic states 
and 3 (dashed lines), and eigenenergies of the undriven po- 
tential (inset), b) Stationary occupations pa in the diabatic 
basis (solid lines) and approximation based on effective rate 
R ac , Eqs. (j6j) and (dotted lines), c) same as b) with loga- 
rithmic axis for pa. The parameters are /1 = 0.08, h c g — 0.1, 
fi = 0.08165, (3 = 60, r) = 10~ 6 and uj c = 100. 
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FIG. 5: a), b) Peak height P r (Ao) of the total probability in 
the right well vs. effective coupling strength F/A and temper- 
ature 1//3. Diamonds indicate the parameters of Fig. [4] The 
dashed line in b) gives the probability in the right well without 
driving, A = 0. The dotted line at 1//3 = E1—E0 indicates the 
transition between the high and the low-temperature regimes. 



and £3 — Eq w 3il in Fig.[4j We have observed switching 
also in the case of strong driving, supporting the gener- 
ality of the proposed switching mechanism. 

IV. CONCLUSION 

In conclusion, we demonstrate a new switching mech- 
anism for an asymmetric double well potential under a 
weak periodic driving and a weak coupling to a heat bath. 
As the origin of the switching we identify an avoided 
crossing in the quasienergy spectrum of the system. Un- 
der its influence the asymptotic occupations of all Flo- 
quet states dramatically change even if both involved 
Floquet states have just small occupations. We explain 
this switching mechanism by an effective rate equation 
at the avoided crossing. 

We now briefly discuss possible advantages of the 
switching mechanism in applications: (i) If one uses a 
laser for the periodic driving, the amplitude dependence 
of the switching mechanism and the beam profile allow 
switching at a 3D spatially localized position with a res- 
olution smaller than the focus width, (ii) In situations 
where a theoretical modeling of the system, e.g. a com- 
plex molecule, is not achievable and no other switching 
mechanism is known, the generic appearance of avoided 
crossings in time-periodically driven systems suggests the 
existence of driving parameters for the desired switching. 

We emphasize that this switching mechanism is com- 
pletely different from standard techniques which allow to 
transfer a wave packet from one well to the other by res- 
onant or near resonant driving and negligible coupling to 
a heat bath. There one has to prepare a specific initial 
wave packet and has to apply the driving for a specific 
duration. In contrast, here the initial state of the system 
is arbitrary, the duration of the driving is arbitrary (if 
larger than the relaxation time) , and the presence of the 
heat bath is essential. 
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